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Oscillations of neutrino emerging from a supernova core are studied. In this extremely high 
density region neutrino self interactions induce collective flavor transitions. When collective 
transitions are decoupled from matter oscillations, as for our chosen matter profile, an ana- 
lytical interpretation of the collective effects is possible, by means of a mechanical analogy 
with a spherical pendulum. For inverted neutrino hierarchy the neutrino propagation can be 
divided in three regimes: synchronization, bipolar oscillations, and spectral split. Our simu- 
lation shows that averaging over neutrino trajectories does not alter the nature of these three 
regimes. 



1 Introduction 



Supernova neutrino oscillations are a very important tool to study astrophysical processes and 
to better understand neutrino properties l^^. When neutrinos leave the surface of the neutri- 
nosphere, they undergo vacuum and matter oscillations. Beside this, in the first few hundred 
kilometers neutrino-neutrino interactions induce collective flavor transitions, whose effect can 
be very important, depending on the neutrino mass hierarchy. Self-interaction effects are ex- 
pected to be non negligible when ;u(r) ~ uj, where ^(r) is the neutrino potential associated 
to the neutrino background {fi = ^/2GF{Ny{r) + Ny{r))^ analogously to the MSW potential 
A = ^/2GFN^-{r) ) and UJ is the vacuum oscillation frequency. We neglect the solar mass square 
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and consider a two-neutrino mixing scenario 



where the oscillations are governed by the mixing angle ^13. Since in the supernova context 
and Uj- cannot be distinguished we generically speak of Ve ^ t^x oscillations. In our work we 
assume Am^ = 10~^ eV^ and sin^ ^13 = 10~^. Figure 1 shows the radial profiles of the matter 
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Figure 1: Radial profiles of the neutrino self-interaction parameter fi{r) — \/2Gf {N + TV) and of the matter- 
interaction parameter X{r) = \PiGf N^- adopted in this work, in the range r £ [10, 200] km. 



potential A(r) and of the neutrino potential /u(r), and the approximate ranges where collective 
flavor transitions of different type occur: synchronization, bipolar oscillation and spectral split. 
The nonlinearity of the self interactions induce neutrino oscillations very different from the or- 
dinary MSW effect. When undergoing collective flavor transition neutrinos and antineutrinos 
of any energy behave similarly, as we will see in the following. This kind of transitions occurs 
for small r, well before the ordinary MSW resonance, allowing for a clear interpretation of the 

numerical simulations. For matter profiles different from our own, the MSW resonance condition 

[2l 

can occur in the same region of the collective transitions: shallow electron density profiles'^ can 
trigger MSW effects around O(IOO) km. In that case it is much more difficult to disentangle 
collective from MSW effects in the results of the simulations. 



2 Reference model and pendulum analogy 

In our work, we use normalized thermal spectra with (Ee) = 10 MeV, (Eg) = 15 MeV, and 
{E^) = (E^) = 24 MeV for i/g, I'ej and Ux, respectively. The geometry of the model, the 
so called "bulb model" has a spherical symmetry, since we assume that neutrinos are half- 
isotropically emitted from the neutrinosphere. Along any radial trajectory there is, therefore, a 
cylindrical symmetry. By virtue of that, we need only two independent variables to describe the 
neutrino propagation and interaction: the distance form the supernova center r, and the angle i? 
between two interacting neutrinos. If the dependence on i) is integrated out, we speak of "single- 
angle" approximation, while the general situation of variable i? is dubbed "multi-angle" case. 
The numerical simulation in the multi-angle case is extremely challenging, since it requires the 
solution of a large system (size of order 10^) of coupled non-linear equations. The propagation 
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Figure 2: Single-angle simulation in inverted hierar- 
chy: Pz (neutrinos) and (antineutrinos) as a func- 
tion of radius, for five energy values. 
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Figure 3: Multi-angle simulation in inverted hierarchy: 
Pz (neutrinos) and P^ (antineutrinos) as a function of 
radius, for five energy values. 



of neutrinos of given energy E is studied through the Liouville equation for the density matrix. 
By expanding the density matrix on the Pauh matrices and on the identity, the equations of 
motion can be expressed in terms of two polarization vectors, P(£^) and P(ii^), for neutrinos 
and antineutrinos, respectively. By introducing a vector B that depends on the mixing angle 
013, and a vector D = J — J that is the difference between the integral over the energy of P and 
P, the equations of motion can be written as 



P = (-FljB + Az ^D) X P , 
P = (-wB -h Az ^D) X P . 



(1) 
(2) 



In the general case, the polarization vectors depend also on the neutrino emission angle (the 
neutrino incidence angle 'd can be expressed in terms of r and of the emission angle at the neu- 
trinosphere ^o)- The electron neutrinos survival probability Pee is a function of the polarization 
vector, Pee = 1/2(1 + Pj/PJ), where the i and / refer to the initial and final state respectively 
(analogously for antineutrinos). The equations of motion for P(-E') and P(-E') can be reduced 
(under reasonable approximations^ to the equations of motion of a gyroscopic pendulum, a 
spherical pendulum of unit length in a constant gravity field, characterized by a point-like mas- 
sive bob spinning around the pendulum axis with constant angular momentum. The pendulum 
inertia is inversely proportional to /i(r), while its total angular momentum depends on the differ- 
ence of the integrated polarization vectors J and Th e motion of a spherical pendulum is, in 
general, a combination of a precession and a nutation'^'^. In the case of normal hierarchy of the 
neutrino mass spectrum the pendulum starts close to the stable, downward position and stays 
close to it, as [x slowly decreases and no collective effect is present. In the inverted hierarchy 
case, the pendulum starts close to the "unstable," upward position. At the beginning, for small 
r, when [i is large (m is small), the bob spin dominates and the pendulum remains precessing in 
the upward position conserving angular momentum^, a situation named synchronization ^ 4!_ 
Nevertheless, since decreases with r, at a certain point any ^13 7^ triggers the fall of the 
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Figure 4: Single-angle simulation in inverted hierar- 
chy: modulus and z-component of J and J. 



Figure 5: Multi-angle simulation in inverted hierar- 
chy: modulus and z-component of J and J. 



pendulum and its subsequent nutations, the so called bipolar oscillations. The increase of the 
pendulum inertia with r reduces the amplitude of the nutations, and bipolar oscillations are ex- 
pected to vanish when self-interaction and vacuum effects are of the same size. At this point, at 
the end of the bipolar regime, self-interaction effects do not completely vanish and the spectral 
split builds up: a "stepwise swap" between the Ue and energy spectra. The neutrino swap- 
ping can be explained by the conservation of the pendulum energy and of the lepton number^^. 
The lepton number conservation is related to the constancy of = Jz — Jz, that is a direct 
consequence of the equation of motion. For a detailed description of the pendulum analogy and 
of our reference model the reader is referred to our previous work^ and references therein. 



3 Simulations 



Figures 2 and 3 show the third component of P and P, as a function of the radius, for different 
energy values, for the single- and multi-angle simulations, respectively. Bipolar oscillations 
starts at the same r and their periods are equal for both v and V at any energy, confirming the 
appearance of a self-induced collective behavior, in the single- and in the multi-angle case. The 
behavior of each Pz and Pz depends on its energy. For neutrinos, Figure 2, the spectral split 
starts around the critical energy Ec — 7 MeV: the curve relative to E < ends up at the same 
initial value {Pee = !)> while the curves for E > Ec show the Pz inversion (Pee = 0). Neutrinos 
with an energy of ~ 19 MeV do not oscilate much, because this is roughly the energy for 
which the initial z^e and i^x fluxes are equal. For antineutrinos, all curves show almost complete 
polarization reversal, with the exception of small energies (of few MeV, not shown in Figure 3). 
Figures m and [5] show the evolution of J and Jz for neutrinos and antineutrinos, in the single- 
and multi-angle cases. The behavior of these vectors can be related to the gyroscopic pendulum 
motion. At the beginning, in the synchronized regime, all the polarization vectors are aligned so 
that J = Jz and J = J z'- the pendulum just spins in the upward position without falling. Around 
~ 70 km the pendulum falls for the first time and nutations appear. The nutation amplitude 
gradually decreases and bipolar oscillations eventually vanish for r ~ 100 km. At the same time, 
the spectral split builds up: antineutrinos tend to completely reverse their polarization, while 
this happens only partially for neutrinos. As said before, also for antineutrinos there is a partial 
swap of the spectra for ~ 4 MeV. From Figure [5] it appears that bipolar oscillations of J and J 
are largely smeared out in the multi-angle case. The bipolar regime starts somewhat later with 
respect to the single-angle case, since neutrino-neutrino interaction angles can be larger than 
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Figure 6: Single-angle simulation in inverted hierarchy: final fluxes (at r = 200 km, in arbitrary units) for different 
neutrino species as a function of energy. Initial fiuxes are shown as dotted lines to guide the eye. 

the (single-angle) average one, leading to stronger self-interaction effects, that force the system 
in synchronized mode slightly longer. However, just as in the single-angle case, the spectral split 
builds up, J z gets finally reversed, while the difference Dz = Jz — J z remains constant. Figures[6] 
and [7] show the final neutrino and antineutrino fluxes, in the single- and multi-angle simulations. 
The neutrinos clearly show the spectral split effect and the corresponding sudden swap of fe 
and Vx fluxes above Ec — 7 MeV. In the right panel of Figure [6l the flnal antineutrino spectra 
are basically completely swapped with respect to the initial ones, except at very low energies, 
where there appears an "antineutrino" spectral split. This phenomenon can be related to the 
loss of J and of \Jz \ • Also in the multi-angle case of Figure [7] , the neutrino spectral swap at 
E > Ec — 7 MeV is rather evident, although less sharp with respect to the single-angle case, 
while the minor feature associated to the "antineutrino spectral split" is largely smeared out. 

4 Conclusions 

We have studied supernova neutrino oscillations in a model where the collective flavor transitions 
(synchronization, bipolar oscillations, and spectral split) are well separated from the MSW 
resonance. We have performed numerical simulations in both single- and multi-angle cases, 
using continuous energy spectra with significant iz-'U and z^e-z^x asymmetry. The results of the 
single-angle simulation can be analytically understood to a large extent by means of a mechanical 
analogy with the spherical pendulum. The main observable effect is the swap of energy spectra, 
for inverted hierarchy, above a critical energy dictated by lepton number conservation. In the 
multi-angle simulation, the details of self-interaction effects change (e.g., the starting point of 
bipolar oscillations and their amplitude), but the spectral swap remains a robust, observable 
feature. In this sense, averaging over neutrino trajectories does not alter the main effect of 
the self interactions. The swapping of neutrino and antineutrino spectra could have an impact 
on r-process nucleosynthesis, on the energy transfer to the shock wave during the supernova 
explosion and on the propagation of the neutrinos through the shock wave. From the point of 
view of neutrino parameters, collective flavor oscillations in supernovae could be instrumental 
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Figure 7: Multi-angle simulation in inverted hierarchy: final fluxes (at r = 200 km, in arbitrary units) for different 
neutrino species as a function of energy. Initial fluxes are shown as dotted lines to guide the eye. 

in identifying the inverse neutrino mass hierarchy, even for very small ^13.^ 
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